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A LOWER BOUND ON THE GROUND STATE ENERGY OF DILUTE BOSE GAS 

JI OON LEE AND JUN YIN 

Abstract. Consider an N-Boson system interacting via a two-body repulsive short-range potential V 
in a three dimensional box A of side length L. We take the limit N, L while keeping the density 

p = N/L^ fixed and small. We prove a new lower bound for its ground state energy per particle 



E(N, A) 



>47rap[l-0(pi''3|iogp|3)]_ 



N 

as p — )■ 0, where a is the scattering length of V 



1. Introduction 

The properties of the Bose gas have been studied by many authors O El HI [9], and since the first 
experimental observation of Bose-Einstein Condensation in 1995 [T1, interest in low temperature Bose 
gases are renewed [13l [Ml [TTl [121 UHl [13 [H] ■ One of the most well-known properties of Bose gas is its 
ground state energy in the dilute limit. In this low density limit, the leading term of the ground state 
energy per particle is Airap, where a is the scattering length of two-body interaction potential and p is 
the density. The upper bound for it was first rigorously proved by Dyson 2 , and the lower bound was 
obtained by Lieb and Yngvason [T3]. It is also proved that the leading term is the same in some cases 
where the interaction potential is not purely non- negative [3T] . 

Lee and Yang ^ i8j first predicted the second order correction to this leading term, which is given by 



47ra/9 



128 , 



(1.1) 



This calculation used pseudopotential method and binary collision expansion method [6] . Another deriva- 
tion was later given by Lieb [3] using a self-consistent closure assumption for the hierarchy of correlation 
functions. 

The upper bound of this second order correction is recently obtained in Yau and Yin's work |20) . 
The matching lower bound, however, has not been proved yet, and in fact, the best result so far is 
inapil - Cpi/i^), which was first derived by Lieb and Yngvason [1^. (See also [15] for details.) It also 
should be mentioned that this second order correction is proved for the high density and weak coupling 
regime recently by Giuliani and Seiringer [?]. 

Lieb and Yngvason [T3] used the "cell method" to find the lower bound for the ground state energy. 
In this method, they first converted the interaction potential into a soft potential with the expense of 
kinetic energy. Unlike Dyson [2], however, they did not sacrifice all the kinetic energy but kept a small 
portion of kinetic energy to apply Temple's inequality . Then, after dividing a large box into smaller 
cubic cells, perturbation theory was applied in each small cubic cell. In this way, the correct leading term 
could be obtained. 

We improve this approach with several new ideas including a new method called 'box doubling method.' 
Some notable advantages of our methods are as follows: 

• We use a combination of perturbation theory and Temple inequality. Also, we are able to preserve 
essentially the full kinetic energy after the replacement of the the singular potential by a smooth 
one. 

• Each small cell can be enlarged so that it contains enough number of particles in average. 
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In this paper we show that the error term becomes —Cp^/'^\ logpj^ with the new method. Our main 
result is stated in Theorem [521 

This paper is organized as fohows: In Section 2, we state our main theorem, introduce key lemmas, 
and prove the theorem. In section 3, we show new type of 'cell method', which is equivalent to proving 
some of the key lemmas. In section 4, we introduce 'box doubling method', which enlarges each cell, 
hence makes it contain more particles in average. Some technical estimates can be found in section 5. 



2. Main result and Key lemmas 

2.1. Definition of the system and main result. Let be a non-negative, smooth, spherically sym- 
metric, and compactly supported potential, satisfying V{x) = if |a;| > i?o for x e R'^. The zero energy 
scattering solution (p satisfies 

-A + ^(x) = 0. (2.1) 

We normalize so that lim|2,|_>.oo ^{x) = 1. The scattering length a of y is defined by 

a := lim \x\ (1 - ip{x)) (2.2) 

as in [O] . 

Hamiltonian of the Bose gas system of N particles in a torus A of side length L is given by 

N N 

HN:=-Y.A,+Y,Vix,-x,). (2.3) 

j=i i<j 

We will consider the thermodynamic limit when N and L approach infinity with the density p — N/L^ 
fixed. Later, we will also let p — > 0. We will use the notation a ^ (3 ii a/ f3 = 0{p'^) ot a/ (3 — 0{\ logpl^"^) 
for some c > 0, and a ~ /3 if a//? — 0(1). For simplicity, we will also use the notation 

x„ = (xi,a;2, • • • ,x„), x = XAr. (2.4) 

Throughout the paper, C denotes a constant independent of p. 

Definition 2.1 (The ground state energy). For given Hamiltonian 

n n 

i7„:=-^A,-f^y(a:.-x,) (2.5) 

in the three dimensional torus A, its ground state energy, E{n, A) is defined to be 

E{n,A) := inf / |V,^(x)pdx + ^ / V{x, - x,mx)fdx]. (2.6) 

The main result of this paper is the following theorem. 

Theorem 2.2. LetV be non-negative, smooth, spherically symmetric, and compactly supported potential, 
satisfying V{x) = if \x\ > Rq, whose scattering length is a. Then, there exists a constant Cq > such 
that 

hm ^i^>Anap{l-C,pi\\ogp\') (2.7) 

as p ^ 0. 
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2.2. Notations. Let eo(K) and (1 — t(k,x)) be the lowest Neumann eigenvalue and eigenfunction of 
(—A + ^V) on the ball of radius k ^ Rq, i.e., 

(-A + ^V{x)){l - tKx)) = eo(K)(l - t(k,x)) (2.8) 
with the boundary conditions 

t{k,x)=0, dT{K,x)=0, if|x|=K. (2.9) 

Note that r is spherically symmetric, since V is spherically symmetric. Some properties of r and cq are 
collected in Section [5l most notably, 

3a 3a ,^ C, ,„ , „, 

-<eo <-(! + -). (2.10) 

Using T, we define a function Wj that will be used for our proof, which will approximate the ground 
state of — Aj + ^{^i ~ ^j) when x is fixed except xj. Let tij be the half of the distance of Xi to 

its nearest particle other than Xj, i.e., 

Uj^t; min \xt - Xk\. (2-11) 
We introduce a particle triple cutoff function Fij (x) at a length scale £o defined by 

F, ,(x) = |^ . (2.12) 

I otherwise 

We also introduce a particle triple cutoff function at a length scale £_i ^ ^o, Gij, which is defined to be 

G. ,(x)J' (2.13) 

I otherwise 

Here, £o and £-i are fixed and satisfy £o ^ £-i ^ Rq. Note that tij ^ tji, and in fact, tij does not 
depend on Xj. The same holds for Fij and G^. 

Let us extend the definition of r so that t{k,x) = if |a:| > k. We define Wj by 

VF,(x) - 1 - 51 [F^J{^}T{£o,x, - Xj) + (1 - i^y(x))G„(x)r(t,,,a;, - x,)]. (2.14) 

Our interpretation of this wave function Wj is as follows: When Xi is the nearest particle of Xj, Wj 
becomes 1 — t{£q, Xi — Xj) or 1 — ritij^xi —Xj), depending on tij. Physically, it corresponds to a situation 
in which the particle Xj is moving and all the other particles are fixed. Any particle Xi other than Xj 
has its own 'effect' on Xj, and the radius of this effect is given by £q or ty, provided that tij > Xj is 
affected only by its nearest particle. 

The definition of Wj shows how we can handle the case when three or more particles are close to each 
other; when a particle Xk is close to Xi and Xj, we do not ignore interaction between Xi and Xj but assume 
that the range of interaction shrinks accordingly. Here, the range of interaction is tij, the size of support 
of T. 

We will use length scales from now on such that 

£^i p-i\\ogp\i, 4 - p"*|logp|~*, (i r-. p~i\\ogp\i+i, £2 ^ p'^log p\~i , e^pi\logp\^. 

(2.15) 

Here, 77 is a small positive number such that < 77 < 1/15. We choose the parameters properly to satisfy 

L^h'£2, £2 = 2^£u (2.16) 

where /i, h' are integers. 

We will prove in Section [5] that there exists a constant Cq < 1 such that 

Wj(x) > 1 -Co > 0. (2.17) 
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Other properties of Wj will be proved in Section [5] as well. 

2.3. Basic strategy of the proof. To prove the main theorem, we use a strategy that consists of fom' 
steps: 

(1) For a given Bosonic wave function ^ g L^(A^), we let $j := W^^'i' for j = 1, 2, • • • , A^. Here, 
$j is well defined, since Wj > 0. Using Wj, we convert the hard potential V into a soft potential 
q with the expense of a portion of local kinetic energy: 

Lemma 2.3. Let 

q{K,x) ■.= eo{K)l{\x\<K) (2.18) 

and 

qij := Fijq{io,Xi - Xj) + (1 - Fij)Gijq{tij,Xi - Xj). (2-19) 
Let X — {xi, X2, • • • , xn)- Then, for any 'ii{x) E L^(A^), 

(*,ffiV*)>V/ \W,\^\Vj'f,{x)\^dx+y2 q,,{x)\'^{x)\^dx, (2.20) 

where ^j{x) := Wf^'i'ix) for j = 1,2, - ■ ■ , N. 

Note that, though qij is the soft potential we want to use throughout the paper, it is not 
symmetric, i.e., qij ^ qji. The kinetic energy term here contains |Wjp|Vj$j(x)p, which is 
associated with an operator Tj = ~W~^V jWjV jW~^ . To use this operator Tj, we need to 
impose boundary condition. Noting that Neumann boundary condition gives the lowest energy 
for a Laplacian operator — Aj, we introduce the following condition on Tj that corresponds to 
the Neumann boundary conditions on — A^: 

Definition 2.4. Let ?7 C A be a box. 

(a) Assume that xi, ■ ■ ■ , a!}, • • • , xn are fixed. We say that tpi^j) satisfies 'Wj-Neumann bound- 
ary condition' in U when il>{xj) € L^{U) and W~^^p{xj) satisfies Neumann boundary con- 
dition in U . 

(b) Let j = 1,2,-- - ,n. Assume that Xn+i,--- ,xn are fixed. We say that 'ip{xn) satis- 
fies Wj-Neumann boundary condition' in U when ^(x„) € L'^{U'^) and, for any fixed 

•^1 ? * * * 7 '^j 5 * * * J "^n 3 

V;j(x,) :=i/;(x„) (2.21) 
satisfies Wj-Ncumann boundary condition in U as defined above. 

When we prove lemmas containing |Tyjp|Vj$j(x)p such as Lemma by using Tj, we will 
only consider cases where ^'^■(x) satisfies VFj-Neumann boundary condition, which can be justified 
as follows: For a given tp{xn) G L'^{U'^), we can find ip{'x.n) G L'^{U"') such that tp satisfies Wj- 
Neumann boundary condition in U for j = 1,2, •• • ,n, and the difference between the kinetic 
energies of ijj and is negligible. To see this, we first change V' into 'tp so that Vj{W~^^{xj)) = 

at the boundary and ip differs from ijj only near the boundary. If this change was small enough 
and made sufficiently near the boundary, then, even though VjV' might have been changed a lot, 
the difference between VjV' and Vj?A could be arbitrarily small after integration. This results in 
the negligible difference between the kinetic energies of ip and ip. 

(2) To use perturbation theory, we divide the torus A into small cubic cells. We need to replace 
the soft potential qij by another soft potential qfj so that particles in different boxes do not 
interact via this soft potential. If we just ignore interactions beyond the boundaries of cubic 
cells, however, it will lower the energy significantly. To resolve this problem, we shift the origin 
of division, u, continuously and take an average of energy with respect to u. More precisely, we 
will consider the following: 
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Let u e [-£i/2,ii/2)^ and let G := £iZ^ Ci A. For A G G, we let A„a be a cubic box of side 
length £i centered at (m + A). Here, u corresponds to the origin of the grid that divides A into 
small boxes A„a. This technique leads us to the following lemma: 

Lemma 2.5. Consider a three dimensional grid with the origin u e A. Assume that this grid 
divides A into small cubic cells of side length £i ^ p^^l logpls^'', < 77 < 1/15 such that 

A - y A„A, Kx, n A,a, = A, ^ A,, (2.22) 

A 

where Au\ is a cubic cell of side length £1 and A is an index for those cubic cells. 

For given u, let dA^x be the boundary of Aux and d{xi,dAux) the distance between Xi and 
dAux in A. Define 




if d{xi,dAux) > 24 for any A 
otherwise 

Then, for any ^/{x) = ^{xi,X2, ■ ■ ■ ,xn) G ^^(A^), 



(2.23) 



Y^f q,^{x)Mx)\'dx>mi{-^rJ2f ^A^)\^i^)\'dx. (2.24) 

(3) In each small cubic cell of side length £1 , we use perturbation theory with half of kinetic energy 
to estimate the energy in the small cell in terms of the number of particles in it. When the small 
cell contains too many particles, however, perturbation loses its validity. In this case, a very 
small portion of energy from the main Hamiltonian will contribute the energy we need and 
it is sufficient in order to estimate the lower bound for the ground state energy. More precisely, 
we can prove the following lemma: 

Lemma 2.6. Let xi,X2, ■ ■ ■ ,xn G A. Let B d A be a box of side length £1. Define E{n,B) as 
the ground state energy of n particle system in B with interaction V , i.e., 

n n 

E{n, B) := inf spec {- XI + H ^(^^ ~ ^i^^ (^'^^^ 
with the Neumann boundary conditions. Let 

f(n)-=/"("-^) */"<2p£? 

' ' \{Ap£l - l)n - {2p£lY ifn> 2p£l ' ^ ' ' 

Let e = p3 I log as in (|2.15p . Fix Xn+i, Xn+2, ■ ■ ■ ,xn outside B and let Xn = (xi, X2, - ■ ■ ,Xn). 
Then, for any ip{xn) G L'^{B^) with (pj :— W~^ip for j = 1, 2, • • • , n, we have 

,E{n,B)ml + ]^j^j^jW,\'\V,(p,{x,,)\''dx^ + {j-^ 



Ana , 



|2 



> f[n){l-Ce)^Ur2. (2.27) 

(4) In the cell method in step (3), the number of pairs in a small cubic cell is not n^/2 but n{n— l)/2. 
This gives an additional factor {1 — p~^£^^)^ which is the ratio between n{n—\) and n^. This error 
factor becomes quite problematic because p~^£^^ ^ | log pj^^^'^'' ^ e. To handle this problem, 
we merge 2^'' adjacent small cubic cells of side length £1 into a larger cubic box of side length £2. 
(Recall that we chose £2 — 2^£i.) For a technical reason, we only merge two cells at a time, and 
we apply perturbation theory to achieve a similar result to Lemma 12.61 for this 'doubled box.' 
When we get to the cell of side length £2, we can obtain the following lemma: 
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Lemma 2.7. Let xi,X2, - ■ ■ , xjy G A. Let A^^ C A &e a cubic box of side length £2 ^ | log pl"^ . 
Divide Ai^ into 2^^ smaller cubic cells of side length ti such that £1 = 2^''^2- Call those small 
cubic cells Bi,B2,-- ■ ,B2ah. Define f as in (j2.26p and let 

~ ^^rt(t-i) ift<pel 

Fix Xn+i, • • ■ ,xjs[ outsidc K^^ and let x^-^ = (xi, X2, ■ ' • , ^n)- Define Af{Bk) as 

n 

N{Bk):^Y.^{x,eBk). (2.29) 
Then, for any ip{xn) e L^(A^'J with 4>j — W^^^tp for j — 1,2, - ■ ■ , n, 

n 2^^ 

|T^,f|V,0,(a;„)pda;„ + ^^ / f{M{Bu)mXr.)\'dx^ 

> fin){l-Cp^\\ogp\)^\ml (2.30) 
Proofs of the key lemmas above are given in Section [3] and Section 21 



2.4. Proof of main theorem. 



Proof of Theorem \2.SX Assmiie that p is sufficiently small. Lemma [2?3l shows that, for any ^ e L'^{K^), 
{^,HM^)>e{^,HN^) + il-e)(Y, |Ty,(x)nv,$,(x)|2dx + ^ / q,j\^{^)fd^) , (2.31) 



where $j = W'^'if and e = p^ \ logpp. 

Recall that in Lemma 12.51 we divided A into small cubic cells A„;^ of side length £i , where u is the 
origin of division and A € G is an index for those cells. From Lemma [2.51 we get 



N 



> e{^,HN^) + il-e)J2 f |M^,(x)P|V,<i>,(x)pdx (2.32) 
+ (l-e)inf(-- Tt)^/ I dxidx2 ■ ■ ■ dxi ■ ■ ■ dxj ■ ■ ■ dxN (/ / qij\'^{:x.)fdxidxj). 



The next step is to divide the kinetic energy into the small cubic cells. Let P be the set of functions 
that assign particles to the small cubic cells, i.e.. 



P ^{p\p: {1,2,- ■■ ,N}^G}. 



(2.33) 
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Let p-\X) = {i : p{i) ^X,l<i<N}. Then, 



^ / |W/,(x)nV,$,(x)|2dx 
N r ^ 

EE/ (ni(^fceA„p(fe))da;,)|iy,(x)nV,<i>,(x)p (2.34) 
EE/ ( n 1(2;*; e A„p(fc))dj;fe ) 



j--PU) = ^ k:p(k)=\ 

Here, the last integral of the right hand side of (|2.34l) represents the kinetic energy in A„a- Likewise, we 
also have 

N 



/ dxidx2 ■ ■ ■ dxi ■ ■ ■ dxj ■ ■ ■ dxNCS^ / qij\^{y^)\'^dxidxj\ 
EE/ f n e A„p(fe))d2;fe) (2.35) 

( /,.-^.„ . . . . E n 



ij':i#J:P(i)=p(i) = A k:p(k} = X 

For a given box i3. Let M{B) be the function that indicates the number of particles in B, i.e., 

A^(B) := ^l(a;, G B). (2.36) 
i=i 

Then, applying Lemma from (1^3^ and we get 

e(vI',iJ^M') + -^ / |W^,(x)nv,<i>,(x)pdx 

+ (l-e)(-7 ^-r)^E / dxidx2---dXi---dXj---dxN(y^^ qrj\-^{x)\'^dxidxj) 

> (l-Ce)^ [ 5]/(AA(A„,))|*(x)pdx, (2.37) 

«1 JA™ ^ 

where / is defined as in (j2.26p . Note that in (j2.37p we only used half of kinetic energy. Combining (|2.3ip , 
([232]) . and (|237)) . we get 



(*,i/A.*)>(i-e)E / r,(x)nV,<i>,(x)pdx+(l-Ce)^ / ^ /(AA(A„,))|vI/(x)pcix. (2.38) 

1 "'A" ^ 

Now we consider larger cubic boxes of side length £2- We let A^g be a cubic box of side length £2, 
which is a union of 2"^'* A^a's, where 6 is an index for those larger cubic boxes. With a shorthand notation 
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X = (xi, CC2, • • • , xn), we obtain from Lemma [2771 that 

^ N „ A r 

|W^,(x)nV,$,(x)prfx+(l-Ce)^ / 5]/(AA(A„,))|vI/(x)pdx 

> il-Ce)^[ ^/(AA(AL,))|vI/(x)pdx, (2.39) 

where / is a convex function such that f{n) = n(n — 1) if n < p^2- So far we have proved that 

(vl/,i/^.*)>inf (1-C6)^ / J2f(^(Kom{^)\'d^. (2.40) 

«2 JA" g 

We are left to minimize J2e Since / is convex and J2gJ^{^'ue) = -^j we can use Jensen's 

inequahty to get, for any xi,X2, - ■ ■ , xn, 

Y^ImKo)) > ^fip4) > ^P^ipel - 1) > (1 - Cpi\logp\')Np£l (2.41) 
e 2 2 

Note that this lower bound is independent of u. Hence 



> (1 - Ce)— 3-iV/9£^ = 47rapiV(l - Ce). (2.42) 



This shows that, when p is small enough, there exists a constant Cq such that 



hm ffllA) >4^„p(i_Coe), (2.43) 

A',L— >oo iV 



which was to be proved. 



□ 



3. Lower bound estimates 

In this section, we prove Lemma 12.31 Lemma 12.51 ^-nd Lemma 12.61 which were used in the proof of 
main theorem. 

3.1. Conversion into a soft potentiaL 

Proof of Lemma\2^ Fix j and consider xi, a;2, • • • , • • • , xat to be fixed. From that $j = W^^'^^ we 
get 



/ |V,4'(x)|2(ix, + [ ^(^» - ^j)\^i^)\'dx, (3.1) 

= I |W^,(x)nv,<i>,(x)pdx, + E / K- W( - + ^^(^« - ^.■))H^.-(x)l |<i>,(x)pdx, 

For each i 7^ j, we have either i^y = 1 or Fij = 0. Note that Fij is independent of Xj. 

(1) When i^y = 1, consider B{xi,£o), a ball of radius £q centered at Xi. When xj £ B{xi,£o)j 
Wj(x) = 1 — T{£o,Xi — Xj), thus 

( - A, + ^V{x, - x,))W,{^) = eo(4)W^,(x). (3.2) 
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When xj ^ B{x^,io), {- Aj + \V{x^ - Xj))Wj{-x.) = 0. Thus, 
W,(^){~ A, + - x,))W,{^)] \^,{^)\'dx, 

eo{io)^\x^ ~ x,\ <4)|W^j(x)n$,(x)|2dx, = / F,,(x)q(4,a;^-x,)|*(x)|2dx,. (3.3) 

JA 

(2) When Fij — 0,Gij = 1, consider B{xi,tij), a ball of radius tij centered at Xi. When xj € 



B{xi,tij), Wi{x) = 1 - T{tij,Xi - Xj), thus 

1 



-A 



-r -V{x, - Xj))Wj{x) = eo{Uj)Wj{x). 



(3.4) 



When Xj i B{x^,UJ), W^j(x)( - A^ + ^V{x^ ~ Xj))Wj{x.) = 0. Thus, 
W,{k){ - A, + V(x, - x,))W,{x)\ \<i>,{x)\'dx, 



(3.5) 



eoiUj)l{\x^ ~ Xj\ < Wr|*j(x)rdT, = / (1 - Fy(x))Gy(x)q(ty,a:, - Xj)\^{x)\'dx^ 



(3) When Fij = Gij — 0, we have that M/j(x) = 1, thus we get 

W^,(x)( - A, + ^Vix, - x,))W,{x)] |<i>,(x)pdx, > 0. 
From cases (l)-(3), p.ip implies 

^ (iv,*(x,)p + i E ^(^'^ - 



(3.6) 



> 



|2|r7 ,f> /'„m2. 



\W,{x)\'\V,^,{x)\'dx, (3.7) 

+ J2 [F^3{^)q{eO,X^ - X j) + (1 - F,j {x))G^j {x)q{Uj , Xi - X j)] 1 4' (x) | ^ . 

To get back to the N particle problem, we first integrate p.7p and summing it over j gives the desired 
lemma. □ 

3.2. Decomposition of A. 

Proof of Lemma\K^ Recah that we let u G [-li/2,li/2f = T,G = l^I? n A, and for AeG, A„a be a 
cubic box of side length l\ centered at (m + A). Here, u corresponds to the origin of the grid that divides 
A into small boxes A„a- Note that the positions of those boxes depend on u. 
Define x«a by 



1 ifxe [-^ + 2£o,- 
otherwise 



(3.8) 



and Xu\ '■= xk^ — w — A). If Xu\[xi) = 1, then Xi G A„a and Xi is not within distance 2£o from the 
boundary of A.a\- Note that 

44,3 



(3.9) 



for any x G A. This means that the probability of having x G A to satisfy XuX = 1 for a A is {£i — A£o)^ /If. 
From the definitions of qfj and Xu\, we have 



9»j(x) = E Xuxix^)qij{x). 



(3.10) 



AeG 
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{qij depends on u, but we omit it.) Thus, for any 4'(x) e L^(A^), 

1^ [ du I g^|^'(x)pdxi(ia;j = / ( du'S^ Xu\{xi)qij\^{y.)\'^dxidxj 
|i I Jr JA2 l-L I Ja2 Jt ^ 

= (^1^)3 j g,^.|vp(x)pdx,dx,. (3.11) 
«1 Ja2 

Integrating p.lip with respect to dxidx2 ■ ■ ■ dxi ■ ■ ■ dxj ■ ■ ■ dx^ and summing over i,j gives 



•^Ja" |i I Jr «i-«o T^JA" 

> M{-^fYl / ^^l*WPrfx. (3.12) 

This proves the lemma. □ 
3.3. Lov^rer bound estimate - Small cubic cell. 



Proof of Lemma 1^.61 We consider the fohowing cases: 
(1) When n < 9p£f : 

When n < 1, f{n) = and the lemma is trivial. Suppose that n > 2. Define an operator 

T, = -Wr^V,WfV,Wr^ (3.13) 

on the functions in L^{B'^) with Wj-Neumann boundary conditions in B in the sense that, for a 
function V' S L^(5"), 

= -l^-iV,(W,2v,(W^rV)). (3.14) 
For i,j G {1, 2, • • • , n}, i 7^ j, we want to consider an operator 

, , . ^1 ^3 



which is defined on functions in L^(i?") satisfying IVi-Neumann boundary conditions and Wj- 
Neumann boundary conditions. 

We first estimate a lower bound for 

, , ^1 ^3 



4(n^l)+(^)^^ (^-^^^ 

by Temple's inequality |19) . with [4(n — \)]^^Tj as the unperturbed part in the first order per- 
turbation theory. 

To find the gap of Tj , we first notice that Wj , which is defined in (|2.14p , is the ground state 
of Tj. For a function Tp{xj) with VFj-Neumann boundary condition in B with 

Wjip{xj)dxj = 0, (3.17) 

B 

we apply Poincare's inequality to obtain 

\ilj{xj)-{l)\'^dxj<C\B\^ I \S/xl}{xj)\^dxj, (3.18) 

JB Jb 

where 

^ = 4t / ^{xj)dx,. (3.19) 



A LOWER BOUND ON THE GROUND STATE ENERGY OF DILUTE BOSE GAS 



11 



We also have 



iV'(xj) - il:\'^dxj = \ijixj)\'^dxj - j ^ i){xj)dxj) 

\4'{x,Tdx, - ^(^^(1 - W,mx,)dx,f > \4,{x,)\''dx, - ( / H{x,)\dx,y 



> (1-c^) / \^Pix,)\'dx,. (3.20) 
Jb 

Thus, we can see that, in a box of side length £i, 

(gap of Tj ) > C^r^- (3-21) 
In order to use Temple's inequality, we first need to check that 



+ (3-22) 



4(n-l) ^£i-Mq 
which is obvious since Tj > and q^j > 0. Let {F)\Yi denotes 

{F)w, = f dx,F\W,\y [ dx,\W^\^ (3.23) 
J B Jb 

We also need to have 

(gapof^)»(^)3(,~).^., (3.24) 
which can be easily seen from that 

(gap of -^^^^) > Cn-H-^ > logpl-t-^" 

» Cpi\logpr'>C\\q,,\\^>ij-^f{^,)w,. (3.25) 
Thus, we can indeed use Temple's inequality to obtain that 

+ ^. _AeJ *J ^ ^. _A,J Wj/w, ~C ——2 . (3.26) 



Hence, 



4(n-l) '£i~Uo' - 'ii-Aio' ' " ' n-H 
We want to estimate {qij)wj £^nd {qij'^)wj ■ It follows from Lemma [5.1l that. for all £-i < k < £o, 

C 

T{K,x,-Xj)<- — 1 (3.27) 

\Xi Xj\ 

-^j I -^j I 

Let S-i and Sq be sets of all points in B that are not within distance 2£^i and 2£o, respectively, 
from any of xi, a;2, • • ■ , x^, • • • ,Xj,--- ,xn, i.e., 

S-i^{xeB:yk^ ij, \x-Xk\ > 2£^i}, (3.29) 

So = {xeB:yk^iJ,\x-Xk\ >2£o}. (3.30) 
Let i? be a set of all points in B that arc not within distance 2£q to the boundary of B, i.e., 

B ^{xe B : d{x, dB) > 24}. (3.31) 
By definition, qfj = if Xi ^ B\S-i or Xi G B\B, and qTj — Qij if Xi G S-i B. 



JI OON LEE AND JUN YIN 

When x^e Son B, 

f ^ C f c 

q^j\Wj\^dxj > / r)dxj>eo{io) (1 - , 

JB Fj— a-jl J\xi-Xj\<io \Xi — Xj\ 

> (1 - ■^)eo(4) / Idx, > Anail ~ (3.32) 

*o J\x,-xj\<eo ^0 

When X, G {S-i\So) n B, 

r Q r (J 

q.^\Wj\^dxj > / - I r)dxj > eo{Uj) / (1 - -)dxj 

JB J\xi-xj\<tij Fi — 

> {l-^)eo{Uj)j ldxj>AT:a{l-^). (3.33) 

f^ij J\xi-Xj\<tij ''ij 

For {qij^)wj , when a;^ € S'o n B, 

^^\Wj\'dxj < [ ^^dx, < Ci^\ (3.34) 
; J B 

and, when Xi £ (5_i\5o) n B, 

j qr/\Wj\^dxj < J ^^dxj < Ct-"^. (3.35) 



Since we know from Lemma 15.21 that 



£?(! - - Cr,') < \W,\^dx, < il (3.36) 



from (p:^ and ([2^, we get 



^ ^J^^/iryi^--)-^i^^^SonB) + {l--)-l{x,£{S^AS^)nB)], (3.37) 
and, from p.34p and p.35p . we get 



)Wj {<li])w, ^ {qij )Wi „2 f ~2iTi7 l2 J if ITI7 l2 7 



B 



3 



< C-gi^r'-U^;, e5onB) + C-3i^r^-i(x, e (5_i\5o)nB). (3.38) 

Inserting p.37p and p.38p into p.26p . we obtain that 



4(n-l) '^i-4£o 

<ir^''^ (,l - £ - cf) ^fe c 5„ n B, + ,1 - I - c|) . !(.. c (5_ ASo) n B,) 

ij^^f^ ( (1 - C^) ■ i{x, e 5o n i?) + (1 - C^) ■ i(x. e (5-i\5o) n B)\ . (3.39) 

Here, for the last inequahty, we used nt\ > t1 and, when x^ G (S'_i\S'o) n S, > > 

/2 > /2 
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Let 

e,(x„) = C^) ■ e ^0 n B) + (1 - C^) ■ e (5_i\5o) nB)\. (3.40) 

Note that Sq and 5-1 are independent of Xi and Xj, and is independent of Xj. We apply 
Temple's inequality to 



4(n-l) '<'i-44 
with [4(n — xy^^^Ti as the unperturbed part. Then, we get 



+ (ir^n^ (3.41) 



We now estimate (^i)vFi- By definition, 

{i^)w.>^l^ I ^^\W,\^dx, 

> Jl'm\'dx.-i[ Ji'^dx.. (3.43) 

*i *o Js-inB *i J(S_i\So)ns Sj 

To estimate the last term in the right hand side, we note that (1) tij = ^\xi — Xk\ for some Xk 
other than Xj and (2) the union of annuli {x : 2£-i < \x — Xk\ < 2io} for all Xk & B other than 
Xi and Xj covers (S'_i\S'o). Thus, 



\x^-Xk\ ^dXr. (3.44) 



<2£„ 



To estimate the first term in the right hand side of (13.431) . we use 

/ ^\W^fdx^> [ ^dx^- [ il-\Wrf)dx^>\S-inB\^Cn£l, (3.45) 
Js-inB Js-inB Jb 

where the last inequality follows from Lemma 15.21 Since IS"-!] > if — Cn{£^i)^ and \B\ = 
{ii - 44)^ from (13^ . (l3:44ll . and (13:45]) . we get 

{ii)Wi 



> 



-1 *1 *■! i^.j._^^J2i-i<\x,-Xk\<2io 



£\ \Xi Xk\ dxi 



>- (^)'^(1-C»<tf -Cf -C„|-C„^U|^,, (3,46, 

Since {^f)w. < \Ml<Ci^^, 

^4Yp2 ^ C'n^r^- (3-47) 
Thus, combining (|33H1), dSHH), ^ZB, and (IgTTT)) . we obtain that 



4(n-l) 4(n-l) '4-44' 
> (l-Cp*|logpP)iy. 
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This implies, when Xn+i, ■ ■ ■ ,xn are outside of B, for any V'(xn) G L'^{B^) with 0; = Wi for 
Z = 1,2, • • • ,n, 

^ ^ (|iy.nv.<^,(x„)p + |VF,f |V,0,(x„)Hrfx„ + (-^i— )3 / 5^|^(x„)pdx„ 



> (3.49) 

This shows how we can get the fower bound for the ground state energy in a soft potential 
regime, qij depends on particles other than Xi and Xj, but overall effect from them is insignificant 
and becomes a small error. 

We apply this argument to all 1 < «, j < ?^, i 7^ j. After summing over all indices i and j, we 
get 

|W^,f|V,</),(x„)pdx„ + (-ii— )3^ / 5^|V'(x„)pdx„ 

> ri{n-l){\-Ce)^\\nl. (3.50) 

(2) When 9/9£? < n < <de-^pll: 

Let p := 9p€f . Here, n satisfies that 

^ 1, np^«l,4«l. (3.51) 



We form particle groups in B, each of which contains at most p particles, such that Gi = 

{xi,X2, ■ ■ ■ ,Xp\, G2 = {Xp+i, ■ ■ ■ ,X2p}, • • • , G^n/p] = {ln/p\-l)p+l i ' ' ' j^J^n/pJ^}, G^n/p]+l = 

{x^n/p\p+i I ■ ■ ■ ,Xn}- Note that the number of groups with p particles is [n/pj and [n/pj > n/2p. 
For i, j e Gfe, we consider 

' +iT^m (3.52) 



4(p-l) 4(p-l) 'h-Uo' 

which is defined on functions in L'^{B") satisfying VFi-Neumann boundary conditions and Wj- 
Neumann boundary conditions. We then use the Temple's inequality as in case (1) to get, for 
any V'(x„) G L^{B'') with 0, = Wf^^j for Z = 1, 2, • • • , n. 



(|M^,nV,0,(x„)p + |W-,f |V,0,(x„)p)dx„ + (^^)=^ / S;|^,(x„)|2dx 



4(p-l)yB„^' " '--J' , w-.v""/, y-"^^ 'ii-Uo 

> (1-C.(^-Gp|-Gn|-Gnpn|^)|i|„||2. (3.53) 

We apply this inequality to all particles in the particle group Gi, then to all particle groups 
in B. Using p. 511) . we obtain 

JE/ |VF,nV,0,(x„)pdx„ + (-^)3f^ / S^|^,(x„)pdx„ 
> L^Jp(p-l)(l-Gnl^-Gp|-G.|-Gnpn!|H)^Mi (3.54) 

We note that, while in Lemma [2.61 the two cases n < 2pi\ and n > 2pif are distinguished, we 
choose the size of p to be 9p^f . If we would choose the size of p to be 2pii, we would get a lower 
bound that is not convex aX n — 2p^2j since the lower bound would increase proportionally to 
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n. We need a space to connect two different formulas so that the resuhing lower bound becomes 
convex, thus the distinction is between n < Qpt\ and n > dpif. 
(3) When n > de'^pil: 

In this case, we only use the term eE{n, B) to prove the lemma and ignore the other terms, 
since they are non-negative. It is known that the ground state energy of n particle system in a 
box of side length £i, 

71 71 1 

E{7i,B)>Ai,a{^)n{l-C{^)^) (3.55) 



when the density {n/li) is sufficiently small and 

n 



^i>C(-3) (3.56) 



(See Theorem 2.4 in [15].) 

In our case, if the box B contains m :— %e^^ pl\ particles, then 

h>Cp-"^\ogp\^ »Cp-T^|logp|^ =C(^)-*, (3.57) 



and the density in this case 



^ =9e-V = 9pi|logp|-3^0 (3.58) 
as p — 0. Thus we can indeed use p. 551) to obtain that 

777 ^ '^Tn 1 

E{m,B)>Ana^{l^C{^)Tr) (3.59) 

To find a lower bound of E(n,B), we form particle groups in _B, each of which contains m 
particles. Since we have [n/mj groups of size m, by superadditivity, 

71 47rn 

Ein, B) > l-\E{7n, B) > (4p£3)„(i _ Ce)-^. (3.60) 

m tl 

Since, 7i falls into one of the above categories, we get the desired lemma. □ 

4. Box DOUBLING METHOD 

In this section, we prove Lemma 12.71 

4.1. Lower bound estimate - Large cubic cell. In order to show Lemma [2.71 we need to prove a 
result analogous to Lemma 12.61 when A^^ , a box of side length £2 is given. We note: 

£i^£2- p-^\ogp\-i (4.1) 

More specifically, for a box B with a side length between £1 and £2, we will show that the energy in B 
with n particles {ti ^ p\B\) retains the form, 

4:na{l~Ce)\B\-'^7i{7i~l). (4.2) 

We can see that the error factor 1/n that comes from the ratio between n(n~ 1) and becomes smaller 
as 71 increases, and eventually it becomes 1/n ^ | logpp <C e when the side length of B becomes £2, 
i.e, p£l 

To demonstrate how to enlarge the size of box while retaining the form ()4.2|) . we first consider a simple 
case where we have only two adjacent boxes Aa and Ab with the same size. Suppose that there are n 
particles in A a U A^. Let ua and ns denote the number of the particles in A^ and A^. We assume the 
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potential energy in and As as UAiriA — 1)/|Aa| and nsins — which has the form in (|4.2p . 

For a e M, we define the Hamiltonian as 

n 

= a^(-A,) + UAiuA - f )/|Aa| + UBins - 1)/\Ab\ (4.3) 

i=l 

When a = 0(no kinetic energy), the ground state energy of this Hamihonian is n{n — 2)/\Aa U Ab\- 
But when a = oo (particles are uniformly distributed in A a U Ab), the ground state energy is equal to 
n{n — 1)/\Aa U A^I, which gives the desired form (|4.2p . 

We will show that instead of — aA^, a small potion of Ti can also guarantee the almost- uniform 
distribution and the desired form (|4.2I) . 

This heuristic argument shows our basic strategy in this section, which we call 'box doubling method.' 
Recall that £2 = 2^ti. In this method, we begin from the first step where we have 2^^ small cubic cells 
of side length £1. We consider (2^''/2) pairs of adjacent boxes, and for each pair that consists of two 
adjacent boxes of same size A^i and A^. As explained above, we can get a lower bound for the energy 
of n particle system in A^ U A^ at expense of small potion of T^'s. In this way, we can effectively make 
the size of each box doubled, since the new 'potential energy term' in A^ U Ab also has the form in 
(|4.2p fwhen the density in A^ U Ab is about p). 

At the end of the first step, or the beginning of the 2nd step, we have 2'^''"^ boxes whose dimensions 
are x £1 x 2£i. In the 2nd step, we consider 2^^~^ pairs of those boxes and perform the above process 
again for all the pairs. Keep using this method. At the beginning of the s-th step, we have 2'^''"'' boxes, 
and after applying the above method to 2'^''^''^^ pairs of boxes, the number of boxes gets halved and 
the size of each box doubled. And the new 'potential energy terms' in new boxes also have the form in 
(|4.2p fwhen the density in new boxes is about p) 

We keep repeating it until the side length of a box becomes ^2^ which is when s = Sh, and wc only 
have one box left. The form of the potential term, (|4.2p remains the same throughout this procedure, 
and it can lead us to the desired result. Lemma [^771 We will make this argument rigorous in this section. 

Before we begin the proof, we introduce definitions that will be used throughout this section. 

• s is a non-negative integer that satisfies 1 < s < 3h, where £2 = 2^£i. We let 

£{s) := 2L^J^i. (4.4) 

i{s) satisfies p^^\ logp|3+'' ^ (.1 < £(s) < £2 ^ P^^\ logp|~^, where < r/ < 1/15. 

This s is a label keeping track of which step we are at. We begin from s = 1 and our method ends 
when s = 3h. 

• Aa and A^ are boxes such that the volume of each box |A^| = \Ab \ = 2'*~^£f and the dimensions 
of A A, Ab, and A a U Ab are either 

(1) Aa = Ab= e{s) X e{s) X e{s), AaUAb= £{s) x £{s) X 2£{s), 

(2) Aa = Ab= £{s) X £{s) X 2£{s), AaUAb= £{s) x 2£{s) x 2£{s), or 

(3) Ka = Ab^ £{3) X 2£{s) X 2£(s), Aa U As = 2£{s) x 2£{s) x 2£{s). 

• M.{A) and M{B) are the functions that indicate how many particles among xi,X2, - ■ • ,Xn are 
in Aa and As, respectively, when Xn+i,Xn+2, ■ ■ ■ ,xn are outside (Aa U Ab), i.e., 

n n 

X(A):=^l(x, gAa), A^(i?) :=^1(x, e As) (4.5) 

1=1 i=l 

Note that M.{A) and M.{B) depend on n though we omitted it. 

Aa are Ab are a pair of boxes at the s-th step. Though we consider only two boxes at a time, note 
that we have 2'^''"''"^ such pairs of boxes in s-th step. 
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Note that the potential term in Lemma \2M i.e., the right hand side of ()2.27|) is ATTa\B\^^ f{n), which 
is different from (|4.2|) . f{n) changes from quadratic to hnear at n ~ 'ipt\- At the s-th step, we use fs 
instead of /, and fs{n) becomes hnear when n > K^, i.e., 

• Define {1 < s < 3h + 1) 



t(t - 1) iit<Ks 
{2K, - l)t - \it> K, 



fs{t)---{7L ... (4-6) 



The definition of fa ensures that fs is continuous and convex. And we define the parameters Kg as 
follows, 

• Let 

:= 2pl\ (4.7) 

and we choose Kg such that 

2Ks - Ks+i » I logp|5 /f?:^ and > 2'-^ pll. (4.8) 

For example, 

= (2 - (1 - -^)| logpl"") • r-'p^l (4.9) 

2 2 

We note that fs with a suitable coefficient is our actual potential energy term in (|4.2I) . When s = 1, 
fs is equal to / in (|2.26p . and fs{t) = t{t — 1) when t < p ■ 2^ pi\, i.e., the density is no more than p. 
(Note that 2^^"^ pi\ is the volume of each box in the s-th step.) 

Our proof of Lemma 12.71 requires the following proposition only, where we consider s to be fixed: 

Proposition 4.1. Let n be an integer and 1 < n < N. Assume that a;„_|_i, x„+2, • ■ ■ , a^iv are fixed 
outside (Aa U A^). Then, for any i;{x„) e L^{{Aa U AbT) with (f>j = V for j = 1,2, ■■■ ,n (Here 

(*^1 ; 1 ' ' ' 5 ^n) ) 7 

|10gp| ^J(A^uAb)" J(AaUAb)" \^A\ 

^ (^-CP~')[ -['^'^f M x^)\'dx„. (4.10) 

J(a^uAb)" |AaUAb| 

Here, the constant C does not depend on s. 

Note that the factor (47ra)~^ can be a general constant C in this proposition and subsequent lemmas. 
We keep this factor, however, in order to use it in the proof of Lemma [2.71 

Proposition 14.11 shows the outcome of box doubling method when it is applied to A^ and As at the 
s-th step. We can prove Lemma 12.71 from this proposition. 



Proof of Lemma \KT[ Recall that we defined h a.s £2 = 2''^i. Recall also that we have a cubic box A^^ 
whose side length is £2 and _Bi, i?2, • • • , B23h are small cubic cells of side length £1 such that IJ^^^ Bk = 
Ai^. We have that x„+i,a;„+2, • • • ,xn are fixed outside A^^ and we let x„ = {xi,X2, • ■ ■ ,x„). We want 
to prove that, for a given ^p{xn) € _L^(A"^) with = W^^^ for j = 1, 2, • • • , n, 

|M^,nV,0,(x„)pdx„-f / /(AA(Bfe))|^(x„)pdx„ 

> J{n){l-Cp^\ogp\)^\\nl (4.11) 
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Here, 



Since 

3.= M2<,l-.,|log.|. (4.13, 



we can keep using Proposition 14. II until we have only one box and its side length is £2- Then, we get 

> {l-Cp^\ogp\)A^a j :^iM|^(x„)|2dx„. (4.14) 



/3' '+i(") l/./„ v2. 

Note that we used an idea similar to (j2.34p in order to convert the sum of kinetic energies over small 
cubic cells into the kinetic energy in a larger cubic cell. 

By definition (|4.8I) . K^h+i > 2^'^p£l = p£l. Together with the definition of / in (|2.26p . it can be easily 
checked that /3/1+1 > /. This proves the desired result, (|2.30p . □ 

4.2. Proof of Proposition l47ll To prove Proposition 14.11 we consider large n case and small n case 
separately. When n is large, we use the following lemma to prove Proposition 14. II 



Lemma 4.2. Suppose that n > K^+i + 2^/ i^s+i • Let iriA, ms = 0, 1, 2, • • • , n. Then, 

mA+mB=n \Aa\ |AaUAb| 

We note this inequality implies ()4.10p directly, since the kinetic energy part is always non-negative. A 
proof of Lemma l4.2l will be given in the next subsection. 



On the other hand, for small n case, i.e., when n < Ks+i + 2y^Ks+i, we are going to prove the 
following inequality: 

|iOgp| ^J(A^uAb)" J(AaUAb)" 1^^-41 

> {l-Cpi)Fs{0,0,n)ml (4.16) 

Here, Fg is a function defined as follows: 

Definition 4.3. For non-negative integers nA,nB, define 

Fs{nA,nB,0) :^ ——— + (4.17) 
|Aa| |Ab| 

To define F(nA, tlb, k), we use the following process: For a fixed fc, consider k particles, j/i, 1/2, • • ■ , j/fc, 
which are uniformly distributed in A^ U A^, i.e., 

Pivr e Aa) - Piy^ e As) = ^ (4.18) 
Let TOyi(fc) and msik) be the number of j/^'s in Ka and A^, respectively, i.e., 

k k 

mA(fe) = ^l(yz G Aa), msik) = ^l{y^ e As) (4.19) 

i=l i=l 
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Extend the definition of F so that 

F,{nA,nB,k) := {Fs{nA + mA{k),nB +mB{k),0))k, (4.20) 

where {■)k denotes expectation with respect to the distribution of t/i, 2/2, ■ ' ■ j (We call those imaginary 
particles yi, y2, - ■ ■ ,yk 'randomized.') 

To prove (|4.16l) . we need the following lemmas that will be proved in subsection 4.3: 



Lemma 4.4. Suppose that n < Kg^i + I^JK^^. Then, 

FMO,n)>il-p)j{s±^ (4.21) 
\Aa U As I 

We note this lemma implies the second inequality of (|4.16l) . For the first one we have: 

Lemma 4.5. Suppose that n < Kg+i + 2^yKs+l■ Fix Xn+i,Xn+2, ■ ■ • ,xn outside {Aa U A^) and let 
Xn ^ {xi,X2, - ■ ■ 1 Xn)- Then, for any ip{xn) G L^((Aa U As)") with 4>j = W~^%p for j = 1, 2, • • • ,n. 



(47ro)-i 
I log p\ 



J2 f \Wj\^\Wj(l>j{xn)\^dx,,+ f Fs{M{A),M{B),0)\^{xn)\''dx„ 



j^-^J{AaUAb)" ^(A^UAij)" 

> il-Cpi)Fs{0,0,n)Ml (4.22) 
Now we are ready to prove Proposition 14.11 



Proof of Proposition \4-l\ When n > Ks+i + 2^yKs+l, the desired result ()4.10p follows from Lemma l4?2l 
When n < Kg+i + 2y/Ks+i, we obtain (|4.16p from Lemma and Lemma which implies (|4.10p . □ 

4.3. Proof of Lemma 14.21 and Lemma 14. 4L 

Proof of Lemma \4-S\ Since fs is convex, 

TL 

min f,{mA) + fsims) > 2/,(-). (4.23) 

mA+mB=n Z 

Thus, it suffices to prove 

2/.(|) > lfs+i{n). (4.24) 



(1) When Ks+i + 2y/Ks+i <n< 2Ks, we have 

4/.(^) - fs+iin) = in' - 2n) ~ [{2X^+1 - l)n - K^^,] 
= {n-Ks+i-\f ~K,+^--^ (4.25) 

With n - Ks+i > 2y/K^, we have it is above zero. 

(2) When n > 2Ks: we compare the derivatives of both sides of (I4.24p . 

Since we also have 

2fsiK,) > i/,+i(2X,), (4.27) 
We can see that, for any n > 2Ks, 

fl 1 

2fs{^) > ^fs+i{n). (4.28) 

Hence, we can get the desired lemma from cases (1) and (2). □ 
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Proof of Lemma \4.4\ Suppose that we have n randomized particles yi,2/2, ■ ■ ' jVn- Each particle is uni- 
formly distributed in Aa U As so that, for i — 1,2, ■ ■ ■ ,n, 

P{y^ e Aa) = P{y. e Ab) = ^ (4.29) 

Let 

n k 

= ^ e Aa), TOi3 e Aij), (4.30) 

1=1 1=1 
and (•) denote expectation with respect to the distribution of j/i, 1/2, ■ • • , i/n- Then, 

Fs(0,0,n) = . (4.31) 

To compute {fs{mA)), we first calculate (m^ — rn^) and estimate the difference {m\ — tua — fs{mA))- 
The former is 

2 2 

{m\ - ruA) - {ruAf + {{m\) - {rriA)^) - (mA) = \ + \-'^^\-\- (4-32) 

To estimate the latter, we use the Chernoff bound (See Corollary 4.9 in [H].) for binomial distribution, 
which becomes, in this case, 

P{mA>n/2 + C)<e-^. (4.33) 
where we used the mean and variance of are n/2 and n/4. Since 



(|4.33l) implies 



{m\ -mA- fsiruA)) = ^("^^ = ' (C - < ^ P{mA > C) ' (C " 



e ^ {C-K,f< e- (C- j^.)^dC < Cn^e-- (4.35) 



To estimate the last term in the inequality above, we use n < Kg+i + 2-^X7+1 and (j4.8p and obtain 



Ks-"^ > i^s-^^- Vi^ITr»|logp|^0^^. (4.36) 



Hence, 



which gives 

Together with (|4.35p we get 
thus, 

Therefore, from that 



„2g-<'^= „"^" < pn^. (4.38) 
(mi-mA-/s(TOA)) (4.39) 



(/s(mA))>(l-p)(---). (4.40) 
ifsirriA) + fsime)) = 2(/.(mA)), (4.41) 
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we obtain, when n < Ks+i + 2yjKs+i, 

^ , , 1 , ^r? ~ n , ^ — n , , f,+i (n) 

FsiO, 0,n) >-,il-p)^ = i^- P)^^^ > (1 - ^'-''^ 

which was to be proved. □ 

4.4. Proof of Lemma 14.51 In this subsection, we consider s to be fixed and let £ = i{s), F = Fg and 
n < Ks+i + 2y/Ks+i. We wiU see that the ground state energy is equal to the expectation with respect 
to the uniform distribution of particles up to a small error. And this is guaranteed by the small portion 
of the kinetic energies ^ Tj . The following lemmas show the idea: 



Lemma 4.6. When n < A's+i + 2y^Ks+i, for any fixed X2, X3, • • • , xn with a;„+i, Xn+2, ■ ■ • , xn outside 
{A A U Ab), assume ua particles among X2, - ■ ■ ,Xn are in A a and ub particles in Ab(Vi^ + ub = n~ 1). 
So, when xi is in Aa the total energy is F(nA + l,nB,0), otherwise it is F{nA,nB + 1,0). Let 

Ti = -VKf ^ViVt^fViWf ^ (4.43) 

defined on functions in L^(A^UAb) with the Wi- Neumann boundary conditions in A^UAs- Then, there 
exists a constant C such that, we have the following operator inequality: 

\, \ Ti + F{nA + 1, ns, 0) • G A^) + ub + 1,0)- G A^) 

> FinA,nB,l)-C'^-^^^^-Cnph-^ (4.44) 

We note that the LHS of (|4.44p depends on xi, but the RHS of (I4.44[) does not. Instead, the function 
F in the right hand side has 1 instead of in its third argument, which means that, when we compute F, 
we need to consider one imaginary particle whose distribution is uniform in A^ U Ab as in the definition 
of F. Thus, we can say that the particle xi got randomized. 

We note that the right hand side has a term —C'{nA — n^)^^"^! logp|. When we apply perturbation 
theory with Xk, 2 < k < n, we also need to take this term into consideration. The following lemma shows 
an outcome of perturbation theory for general k: 

Lemma 4.7. With the assumption as above, except that ua particles among Xk+i, ■ ■ ■ ,Xn are in Aa and 
Ub particles in Ab (riA + ub = ri ^ k), we have 

(47ra)-i . fr.r . . , in ^ui ,AnA + l~riB^^ 



F{nA,nB + 1, fc - 1) - C'ik - 1) ^''^4^Qg''^Pi^^' ) -lixkeAB) (4.45) 
> F{nA,nB, fc) - C'k ^J^f, ~ "f r - Cripk-'. 

^4|l0gp|-l 

Here, C does not depend on k and Tj. is defined as 

Tk = -W^^VkW^VkW^^ (4.46) 
We can prove Lemma 14.51 using the above lemma n times. 
Proof of Lemma Let 



X«(v4)= gA^), xW(S)= 1(2;. e As). (4.47) 

i=k+l i=k+l 

Note that M{A) = M'^°\A) and M{B) = M^°^{B). 



22 



JI OON LEE AND JUN YIN 



From Lemma I4.6[ we get 



(47ra)-i 
I log p\ 



|W^inVi(^i(x„)|2dx„+ / F{M{A),M{B),Q)\^{^n)?d^^ (4.48) 

(AaUAb)" J(AaUAb)" 



J(AaUAb)" * I ^OgPl 



As a next step, we apply Lemma 14.71 with the right hand side of the above equation. Then we get, 
(47ra)-i 



I log p\ 



iw^2riv2</)2(x„)rdx, 

(A^UAi,)" 



> 



[F[M^^\A)M'\BU)-C'^^^^^^^^^^^ (4.49) 
(a^uAb)" ^^logpl 

(a^uAb)" « Mogpl 

Thus, 

E TT^ / |M^.nV.0.(x„)|2dx„+ / F(X(A),X(B),0)|V'(x„)pdx„ 

4=1,2 l^°g^l "'(AaUAb)" J(A^UAb)" 

> / [F(A^(^)(A),A^(^)(B),2)-C'^:^:1^^^^^^^ 

We keep applying Lemma [4.71 Since Ai'^'^^^A) = M.^''^{B) = 0, when k = n, we get 

"(47raK/' |W^,|2|v,<^,(x„)|2dx„+ / F(X(A),A^(i3),0)|^(x„)pdx„ 

|iOgp| J(AaUAb)" J(AaUAb)" 

> / [F{M^-HA),M^-\B),n)-C'^-^^^-^^^^ 
J(a^uAb)" ^^logpl 

= [F{0,0,n)-Cn^ph-^]\\m. (4.50) 

Since f is convex, from Jensen's inequality, 

F{0,0,n)>F{^,^,0)>Cn^r\ (4.51) 

Thus, we have 

TT^E / |M^,f |Vj0,(x„)|'rfx„+ / F(A^(A),A^(i?),0)|V(a;„)pdx„ 

|iOgp| P^J(AaUAb)" i(AAUAB)" 

> (l-Cp3)F(0,0,n)||^(a:„)||2dx„, (4.52) 
which was to be proved. □ 

Since Lemma 14.61 is a special case of Lemma 14.71 we only prove Lemma 14.71 

Proof of Lemma \4- 7\ Let 

^^1^^^^ _ iF{nA + l,nB,k - 1) - F{nA,nB,k - 1) if x^. G (4:53) 
\F{nA,nB + l,fc- 1) - F{nA,nB,k - 1) if Xk G A^ 



and 



AMx.)= il^'^'f. (4.54) 

I [UA - Ub - I) - (UA - ub) it G As 
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Consider 

-Tj rTk+9k-C {k-1)— — - 4.55 

|logp| ^llogpl 1 

with (47ra| log p\)^^Tk as the unperturbed part. We want to use Temple's inequality to get a lower bound, 
i.e., 

^ (4.56) 
|logp| i*\^ogp\ 1 

C 



> F{nA,nB,k) - F{nA,nB,k - 1) - 



--r{nA~nBy -Cnp^e ^ 

|logp|-i 



which implies Lemma 14.71 So to prove Lemma 14.71 it only remains to prove (j4.56p . 

In order to use Temple's inequality on (I4.55p . we first check if the perturbation part is non-negative. 
This is trivial when A; = 1, so we assume that k > 2. Without loss of generality, assume Xk & A^. From 
definition. 



gkixk) = F{nA + 1, riB, fc - 1) - F{nA, ub, fc - 1) 

{F{nA + 1 + mA{k - 1), hb + mB{k - 1), 0) - F{nA + m^(fc - 1), hb + mB(fc - 1), 0))fc_i 
ifs{nA + 1 + mA{k - 1)) fs{nA + mA{k - 1)) 



/fe-i 



(4.57) 



|Aa| |Aa| 
It can be easily checked that 

fsiuA + 1 + mA{k - 1)) - fsiuA + mA{k - 1)) > C{nA + mA{k - 1)). (4.58) 

Thus, 

gkixk) > Cr^iuA + rriAik - l))k-i > Ckr\ (4.59) 
We can prove similarly that gk{xk) > Ckt^^ when Xk G A^. Since n < Cp£^ ^ £\ logpj"^ and Mk < Cn, 



we have 



g, > Ckl-^ » C J''' > (4.60) 
which shows that the perturbation part of (I4.55P is non- negative, i.e.. 

We also know the gap of (47ra| logpD^^Tfc is much larger than the expectation value of gk in the ground 
state Wk, since 

r^\\ogp\-^ » Cn£-^ > \\gk\\^ > {gk)w,- (4.62) 
Hence, we can apply Temple's inequality on (j4.55l) and obtain 

(47ra)-i Mk 
-Tk+ gk~ C (fc - 1) 



logpi ^ 'eVogp\- 

^4|logp|-i^^'^' ^-2|logp|-i 



> {9k)w, -C (fc~ 1)( „4M„„ -C — (4.63) 



£8 1 log pI -2 

where {gkjw^. denotes 



{9k)w, = ( / gk{xk)\Wk\''dxk)/{ / iT^fcpdxfc) (4.64) 



and other expectations are defined similarly. 
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We want to estimate terms in the right hand side of (|4.63|) . In each estimate, we want to compare the 
expectation {■)wk with an expectation with respect to a uniform distribution. Let {gk)i^ denote 

{9k)u - ( / 9kixk) ■ 1 dxk)/{ j 1 dxk) (4.65) 

and other expectations are denoted similarly. It can be easily checked from the definitions that 

F{nA,nB,k - 1) + {gk)ik = F{nA,nB,k). (4.66) 
We have the following estimates for the terms in the right hand side of (|4.63p . 

(1) Using Lemma [Ol we get 

\{9k)w,-{9k)iA<Cnpilr\ (4.67) 
Together with (|4.66p . we have 

{9k)w, > F{nA,nB,k) - F{nA,nB,k - I) - Cnph-"^ . (4.68) 

(2) Using Lemma [Ql we get 

'l<»r)». - ^J^M ^ « (4.6.) 

On the other hand, we simply follow the definitions to see 

= 1- (4.70) 

Hence, we get 

(3) Using Lemma [Ol we get 

\{9l)w,-{9l)iA<Cn^ptlr\ (4.72) 

\{9k?w,-{9k)l,\<Cn^pilr^. (4.73) 
Thus, since n < Cpi^ < £\ logp|"\ 

To estimate (3^)1^ — {9k} 1^, we again follow the definition and get 

{9l)u - (3fe>L = \ [FiuA + 1, ns, fc - 1) - F{nA, ns + 1, fc - 1)]'. (4.75) 

Suppose that we have (k — 1) particles randomized, namely yi, j/2, ■ • • , J/fe-i as in the definition 
of i^, Km - liTMi . Then, 

[F{nA + 1, UB, fc - 1) - F{nA, + 1, fc - 1)] ^ (4.76) 



(4.71) 



= [{F{nA + 1 + mA{k -l),nB+mB{k- 1), 0) - F{nA + mA{k - 1), + 1 + mB(fc - 1), 0))fe 
< ( [F{nA + 1 + m^(fc — 1), ub + mB{k — 1), 0) — F{nA + ruA (fc-l),nB + l + mi3(fc-l),0)] )^,_^ 
It can be easily checked that, for any non- negative integers mi and m2. 



2 ^ 4(toi - 7712) 



[F{mi + 1, m2, 0) - F{mi,m2 + 1, 0)]^ < ^ \, , (4.77) 
Thus, 

[F{nA + 1, ns, fc - 1) - F{nA, + 1, fc - 1)] ^ 

< A\kA\-^{[{nA+niA{k-l))-{nB+mB{k-l)f])^_^. (4.78) 
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From the definition, 

( [{UA + m^(fc - 1)) - {ns +mB{k~ I)?]) = (ua - ub? + {k-l). (4.79) 
Thus, from (jiJij) . (|4J5l) . gjE]), and (|4J9l) . we obtain 

^-2|logp|-i £-2|iog^|-i ^^^4|iog^|-i +^"^^^ • (4-80) 

(4) Using Lemma [5751 from that n < Cp£^ ^ ^| logp|^^ we get 



u2 



- (M,)^. (M|)i, - (A4)i^ n'pil 



^ .611 " 1-3 « "-"^o^" • (4-81) 



£6 1 log p| -3 ^6|logp|-3 - £6 1 log ^1 -3 

From the definition of M^,, it can be easily checked that 

(M|)i, - {Mk)l^ = iiuA - UB?. (4.82) 

Thus, 

£6|logp|-3 £6|log^|-3 -^6|logp|-3 V ^ 

Here, the last inequality follows from 

n^r^llogpl^ < Cp^£^|logp|3 <pi (4.84) 

Hence, 

f6|logp|-3 

Applying estimates (|4768l) . (jiTTTI) . (|4780l) . and (|4785l) to (|4763l) . we obtain (jiTSe)) . This proves the desired 
lemma. □ 

5. Properties of two body problem and the approximation to the ground state 

Lemma 5.1. Let cq and (p he the lowest Neumann eigenvalue and eigenfunction on the ball of radius k, 
i.e., 

(_A + 1^)0 = 600 (5.1) 

with the boundary condition 

(/)(a;) = 1 and d<j){x) = i/ = n. (5.2) 
Suppose that V{x) = when \x\ > Rq. Then, if Rq, there exist constants cq and c\ such that 

Co < -/-(a;) < 1, l-0(a:)<^. (5.3) 

\x\ 

Moreover, we have 

3a 3a ,^ C, 

-<eo <-(! + -). (5.4) 
Proof. For a proof of the first part of the lemma, 

co<^ix)<l, l-cj,{x)<^, (5.5) 

\x\ 

see Lemma A.l in |5|. To prove the second part, we extend the definition of (j) so that (p{x) = 1 if \x\ > k. 
It is well known that 

(t){x) ( -A + \v{x)) (l){x)dx > Alia. (5.6) 



26 



JI OON LEE AND JUN YIN 



Thus, < 1 implies that 

eo>4na([ \(j){x)fdxy^>^. (5.7) 

Upper bound for eg is also proved in Lemma A.l in [^. □ 

Lemma 5.2. Define Wj as in (|2.14p and x — (xi, X2, ■ • • , xat). Suppose that Xj is in Ai, a box of side 
length £ and xi, - ■ ■ , x} , • • • ,xn are fixed. 

(1) There exists a constant Cq < 1 such that 

Wj(a;) > 1 -Co. (5.8) 

(2) Let £ ^ £q. If Ag contains n particles including Xj, then 

£^{l - Cn-^ - C£-^) < / \Wj\^dxj < (5.9) 

Proof. From Lemma l5.1i we can see that there exists a constant cq that does not depend on k such that 
T^KjXi — Xj) < Co whenever £_i < k < £q. It is clear that 

F,j{x)T{£o,Xi - Xj) > or (1 - F^j{x))Gij{x.)T(tij,Xi - x,) > (5.10) 

implies that Xi is the nearest particle of Xj. Thus, there exists an index k such that 

Wj{x) = 1 - [Ffej(x)r(4, Xk - Xj) + (1 - Fkj{x))Gkj{x)T{tkj,Xk - Xj)], (5.11) 

hence, we obtain the first part of the lemma. 

To prove the second part of the lemma, we note that Wj = 1 unless there exists a particle Xk such that 
\xk — Xj\ < min{£o, tkj}. Assume that there exists such k. Consider first a case where Xk £ Ag. When 

tkj > £o, 

Wj\^dxj>[ (l-| :fdxj> I dxj-C£l, (5.12) 

\xk-Xj\<eo J\xk-Xj\<eo \Xk — Xj\ J\xk~Xj\<ea 

and, when t^j < £o, 

mi'dx, > [ (1 - -J—-)'dxj > [ dxj - Gil > --Ktl - C£l (5.13) 

\xk-Xj\<tij J\xk-Xj\<tij Xj\ J\xk-Xj\<tij 'J 

Now consider the other case where Xk ^ A^. Note that dixk.Ai) < £q in this case, where d{xk,Af) 
denotes the distance from Xk to the box A^. When tkj > £o, let Bk be a ball of radius £o centered at Xk- 
Bk n dAg, the intersection of Bk and the boundary of A^, is a disk, and we let be the radius of the 
disk. Since 

If 1 

-dxj < / -j -dxj = 2Trrl, (5.14) 



B^nAi - Xk\ J^/il-rl<\xj-x^\<eo Fj" ~ 

we can see that 

/ mfdxj> [ (1- I ^ y dx,> [ dx.-Crl (5.15) 

JBknAi JB^nAi \Xk Xj\ JBknAi 

We can get the same estimate when tkj < ^o, in this case, Bk becomes a ball of radius tkj centered at 
Xk. Note that all these Bk's are disjoint. If we sum over all such fc, we get 



E / m\'dx,> dx,-G "-i (5.16) 

k:Xkt^t,d{^k.Ai)<eo •'Bkr\Ai k:XktAe,d(xkM)<io k:XkiAt,d(xk.At)<to 

Since Bk n dAg lies on dAg, ^ r\ cannot exceed the area of 9A^, thus. 



E / Wj?dx,> J2 / dx,-G£\ (5.17) 

k:Xk^Ae.dixk.Ae)<eo Bk^At k:x^!f,Ai,d(xk-Ai)<to ^"^^^ 
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Altogether, we get 

[ \Wj\'^dxj> [ dxj-Cntl-Cf = l^{l-cJ^~Cr^). (5.18) 

J Ai JAi ^ 

Finally, Wj < 1 implies that 

/ \Wj\^dXj <f. (5.19) 

JAf 

This proves the second part of the lemma. □ 

Lemma 5.3. Let Aa and Ab be as in Section^ Suppose Xi,X2,- ■■ ,Xn & Aa U Ab and n — 0(p|A^|). 
Let 1 < j < n and let Wj be defined in (j2.14p . Fix xi,X2, ■ ■ ■ ,Xj,--- , xn with Xn+i, Xn+2, • • ■ , outside 
{A A U Ab)- For a given function h{xj) with 

hA if x-i G A A 

h{xi) = / , ' A . (5-20) 
yhB ij Xj e Ab 

let 

{h)w, - / h{x,)\W,\'dx,/ [ W.^'dx,, (5.21) 

{h)i^^ h{xj)-ldxj/ Idxj. (5.22) 

JAaUAb JAaUAb 



Then, 
Proof. Let 



Then, 



\{h)w, - {h)ij < Cpilina^{hA,hB}. (5.23) 



PA^{ \W,\'dx,)/{ m'dx,), (5.24) 

"'Aa JAaHAb 

PB = i[ m\^dxj)/{[ mi^dxj). (5.25) 

JAb JAaUAb 

{h)wj =PAhA+PBhB- (5.26) 

1 nf'^ ^ 

PA^^ + O(^) + 0(^1) = - + 0{p£l) + 0{r'), (5.27) 
and the same estimate holds for ps- Since 

{h)u=lhA + hB, (5.28) 

and <C pi^, we get 

\{h)w, - (Ml J < CpilhA + CpilhB < Cpilm£ix{hA,hB}, (5.29) 

which was to be proved. □ 



From Lemma [5 
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